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We propose a spectral quark model which can be applied to low energy hadronic physics. The 
approach is based on a generalization of the Lehmann representation of the quark propagator. We 
work at the one-quark-loop level. Electromagnetic and chiral invariance are ensured with help of 
the gauge technique which provides particular solutions to the Ward-Takahashi identities. General 
conditions on the quark spectral function follow from natural physical requirements. In particular, 
the function is normalized, its all positive moments must vanish, while the physical observables 
depend on negative moments and the so-called log-moments. As a consequence, the model is made 
finite, dispersion relations hold, chiral anomalies are preserved, and the twist expansion is free 
from logarithmic scaling violations, as requested of a low-energy model. We study a variety of 
processes and show that the framework is very simple and practical. Finally, incorporating the 
idea of vector-meson dominance, we present an explicit construction of the quark spectral function 
which satisfies all the requirements. The corresponding momentum representation of the resulting 
quark propagator exhibits only cuts on the physical axis, with no poles present anywhere in the 
complex momentum space. The momentum-dependent quark mass compares very well to recent 
lattice calculations. A large number of predictions and relations, valid at the low-energy scale of the 
model, can be deduced from our approach for such quantities as the pion light-cone wave function, 
non-local quark condensate, pion transition form factor, pion valence parton distribution function, 
etc. These quantities, obtained at a low-energy scale of the model, have correct properties, as 
requested by symmetries and anomalies. They also have pure twist expansion, free of logarithmic 
corrections, as requested by the QCD factorization property. 
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I. INTRODUCTION 



The term chiral quark model has become a generic 
name for any relativistic field theory aiming at the de- 
scription of the non-perturbative features of QCD. Nu- 
merous approaches (for reviews see, e.g., 0, 0, 0, 0, 
[S Hi and references therein) share a number of 
common features. Firstly, they incorporate dynamical 
quarks as the only explicit degrees of freedom. Secondly, 
they provide particular solutions to the chiral and elec- 
tromagnetic Ward-Takahashi identities. Although there 
is no doubt that chiral quark models provide a reason- 
ably accurate quantitative description of hadronic prop- 
erties, there is a lack of systematics in the construction of 
any particular dynamical model. A source of ambiguity 
is the fact that chiral quark models are, supposedly, an 
approximation to the low-energy non-perturbative QCD 
dynamics. In this regard, an essential ingredient is the 
introduction of a practical suppression of high-energy de- 
grees of freedom, necessary in order to separate the low- 
energy regime, where the model is supposed to work, and 
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the high-energy regime, where the genuine QCD dynam- 
ics, in terms of explicit quarks and gluons, should set in. 
This defines a certain scale, or cut-off, which acquires a 
physical meaning and which should be kept throughout 
the calculation. The precise way how this high-energy 
cut-off should be introduced is not at all clear, although 
some constraints can be imposed on the basis of the rela- 
tivistic and gauge invariance. For a discussions of various 
popular approaches and associated problems see, e.g., 0. 

In this work we introduce a novel approach, the spectral 
regularization of the chiral quark model, based on the 
formal introduction of the Lehmann representation [t| 
for the quark propagator, 

S(p) = I toM-, (1.1) 

where p(u>) is the spectral function and C denotes a con- 
tour in the complex u> plane chosen in a suitable way. 
As will become clear in Sect. IIIII the spectral regu- 
larization allows to explicitly solve the chiral and elec- 
tromagnetic Ward-Takahashi identities in a rather sim- 
ple manner, through the use of the so-called gauge tech- 
nique [l(J, UJ$ ■ The method has already been sketched in 
a previous work by one of us |l2j| . Of course, any solu- 
tion of the Ward-Takahashi identities cannot be complete 
due to the existence of transverse terms, which necessar- 
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ily appear in the underlying theory, and which can only 
be uniquely determined in QCD. In fact, any possible re- 
alization of a chiral quark model, if properly regularized, 
represents in a sense a particular solution to the Ward- 
Takahashi identities. This point is not fully appreciated. 
In Sect. 1 1 1 1 1 we will provide a minimal set of solutions 
to the Ward-Takahashi identities and study their conse- 
quences. 

Throughout the paper the model is considered at the 
one-quark-loop level and in the chiral limit of the vanish- 
ing current quark masses. Therefore our present predic- 
tions are made at the leading- N c level and in the chiral 
limit. 

The proper normalization and the conditions of finite- 
ness of hadronic observables are achieved by requesting 
what we call the spectral conditions for the moments of 
the quark spectral function, p(uj), namely 

Po = J dujp(w) = 1, (1.2) 
p n = J dwcj n p(u) = 0, for n= 1,2,3, ... (1.3) 

It will be shown that the physical observables are pro- 
portional to the inverse moments, 

p. k = J diuuj~ k p(uj), for k = 1,2,3,... (1.4) 

as well as to the "log moments", 

P ' n = J ^lo^/^wVM 

= J dui\og(uj 2 )uj n p(Lu), for n = 2,3,4,... (1.5) 

Note that the conditions i|1.3fl remove the dependence on 
the scale p in l|1.5|l . thus we can drop it. No standard 
requirement of positivity for the spectral strength, p(to), 
is made (see Sect. IIID|l . 

In the present work we do not intend to determine 
the quark spectral function from "first principles" , but 
rather look for general consequences and implicit rela- 
tions which follow from the approach. 

The model with the spectral regularization 
(|l.lll.2ll.3fl . supplied with couplings obtained via 
the gauge technique, possesses simultaneously the 
following features: 

1. Gives finite values for hadronic observables, which 
can be used to fix the inverse moments, (|1.4|) . and 
the log moments, <|1.5ll . 

2. Satisfies by construction the electromagnetic and 
chiral Ward-Takahashi identities, thus reproducing 
all the necessary symmetry requirements. 

3. Satisfies the anomaly conditions. 



4. Complies to the QCD factorization property, in the 
sense that the expansion of a correlator in a large- 
momentum Q is a pure twist-expansion involving 
only the inverse powers of Q 2 , without the log Q 2 
corrections. Thus, the model can be used to com- 
pute the soft matrix elements needed to describe 
the deep inelastic inclusive and exclusive processes 
and to analyze the pion structure function, pion 
distribution amplitude, etc. 

The fact that all above features can be satisfied simulta- 
neously in a chiral quark model is far from trivial . 

One of the advantages of our regularization is that all 
calculations can be directly undertaken in the Minkowski 
space, although nothing prevents us from working in 
the Euclidean space. When dealing with bound-state 
problems, the continuum Euclidean formulation encoun- 
ters practical difficulties in momentum space since going 
to the bound state pole requires a continuation to the 
Minkowski region, and hence a continuation of the quark 
propagator to the complex plane is required. This is a 
problem since much of our phenomenological insight is 
based on the behavior of the quark propagator in the 
Euclidean region. In our approach the analytic contin- 
uation from the Euclidean to the Minkowski space can 
be done in a straightforward manner because the gen- 
eralized Lehmann representation l|l.l|l implies a definite 
analytic structure. 

In the language of the practitioners of chiral quark 
models such as the Nambu-Jona-Lasinio model, the reg- 
ularization introduced by the gauge technique is very 
special because not only it makes the theory finite, but 
also corresponds to taking the infinite cut-off limit in 
those observables which do not depend on the constituent 
quark mass. This includes the proper fulfillment of the 
anomalies. As we will see below, this is a very rewarding 
aspect of the present investigation, which avoids the ar- 
tificial separation between the real and imaginary parts 
in the Euclidean action, or, equivalently, the normal and 
abnormal parity processes in the Minkowski space. At 
the same time, the study of several processes in the high- 
energy limit turns out to be compatible with factorization 
of amplitudes into hard and soft pieces in the twist expan- 
sion precisely because of the regularization and the set of 
conditions (|1.2ll.3fl imposed upon the spectral function 
p(uj). Although such a behavior is expected in pertur- 
bative QCD, it is very difficult to comply to it in tra- 
ditional chiral quark models. This point has been re- 
cently discussed by one of us in Ref. Q, where it is 
pointed out that for the process 7* — ► 7r°7, which in- 
volves the transition form factor F 7 « 7r o 7 (Q 2 ), there is a 
conflict between the chiral anomaly normalization condi- 
tion for F 7 » w o 7 (0) = l/(47r /tt), and the expected QCD 
factorization at large momenta, Q 2 F J , 7r °-y(Q 2 ) — * 2/ w . 
The conflict persists in any standard approach; to fulfill 
the anomaly the absence of a regulator is required, but 
to achieve factorization an explicit regularization must be 
considered. Our model is free of such contradictions, and 
both conditions turn out to be satisfied simultaneously. 
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Actually, as we show in Sects. IVIIII and Ixl for the pion 
the leading twist contribution to the parton distribution 
function (PDF), V w (x), and the parton distribution am- 
plitude (PDA), <Pn{x) the following remarkable relation 
holds at the model's scale Qo- 



(p n (x,Qo) = ^V„(x,Q ) = 1. 



(1.6) 



We stress that the interpretation of results of low- 
energy models in the context of the high-energy processes 
necessarily involves QCD evolution to account for loga- 
rithmic pcrturbative radiative corrections. The soft ma- 
trix elements computed in the model are obtained at the 
working energy scale of the model, Qq, typically quite low, 
and have to be evolved to the scales realized experimen- 
tally with help of the QCD evolution. Only then the com- 
parison to data can be made |13lll4| . In that sense, chiral 
quark models provide the initial conditions for the QCD 
evolution. The phenomenological success of Eq. 11.6fl af- 
ter evolution has been described elsewhere |^. ll5lll6lll7j . 

In Sect. IXIII we provide a model for the moments (|1.4|) 
of the spectral function based on the vector-meson dom- 
inance (VMD) in the pion electromagnetic form factor. 
Actually, only the even and negative moments are de- 
termined by such a method. Remarkably, the positive 
moments obtained by analytic continuation automati- 
cally fulfill the spectral conditions, Eq. (|1.2I1.3|) . and the 
log-moments (|1.5|1 can also be determined. Interestingly, 
the inverse-moment problem for the VMD-inspired model 
can be solved, yielding a simple function with a certain 
cut-structure in the complex w-plane. As a result, the 
model becomes fully explicit and some further results 
can be obtained. Using the insight provided by the VMD 
model we work out in Sect. IXIIII the quark propagator, 
which possesses a certain cut structure but has no poles 
in the whole complex plane. The quark mass function 
agrees remarkably well with recent lattice data 0, 0] . 
We illustrate the power of the method in Sect. IXIVl by 
presenting some further predictions based on the VMD 
model, namely, the pion transition form factor, the pion 
light cone wave function, the non-local quark condensate, 
and the unintegrated parton distribution function of the 
pion. 



II. QUARK PROPAGATOR 

A. Generalized spectral representation for the 
quark propagator 

Our starting point is the definition of the quark prop- 
agator. In the momentum space we have 

S(p) =-i [ d^xe^iOlTUixMO)} |0). (2.1) 



We assume the spectral representation <|1.1[) for the quark 
propagator, where p(uj) is the spectral function, and C 



is a contour in the complex w— plane chosen in a suitable 
way. We do not specify explicitly what the contour C is, 
hence the representation l|l.ll) is a generalization of the 
standard Lehmann representation [9(. To see the con- 
nection, let us consider the special example of a contour 
C running from — oo under the real negative axis and 
crossing through zero above the real positive axis going 
to +oo, yielding the form proposed in Ref. ^} 



S( P ) 



p(uj) 



j> — u) + ie(to) ' 



(2.2) 



where e(uj) = + sgn(o;). Squaring the denominator in 
Eq. (|2.2|) yields a more customary form of the Lehmann 
representation [jj, 



s{p) = [ *M. 

Jo V 



t>pv(u) +ups{uj) 



i0+ 



(2.3) 



where the vector and scalar spectral functions 



p v (u) 



p(uj) + p(-uj), 



(2.4) 



arc defined for ui > respectively. Notice that the func- 
tions pv(ui) and ps(w) are independent of each other. 
Equivalently, the values of p(ui) for positive and negative 
values of u> are also independent of each other. The quark 
propagator in the standard Lehmann representation has 
an analytic structure of poles and cuts on the real axis 
of the complex p 2 — plane, where the positivity conditions 
Pvi^) > and pv(v) > ps{^) for w > hold. With 
Eqs. I|2.4I) . the positivity conditions are equivalent to say- 
ing that p(u>) > for any u>, positive or negative. They 
follow from the requirement of a physical Hilbert. We will 
see in Sect. Ill Dl that chiral symmetry breaking, together 
with the finiteness of hadronic observables, implies that 
in our case a real p(u>) cannot be positive definite, which 
is a simple consequence of the conditions (|1.2ll.3f) if the 
contour is taken to be the standard one. Actually, the 
particular realization proposed below in Sect. IXTJ based 
on the vector-meson dominance of the electromagnetic 
pion form factor shows that one needs in fact a non-trivial 
contour to avoid end-point singularities. We refrain from 
speculating on the possible connection of this fact to an 
indication of quark confinement, nevertheless it is cer- 
tainly true that a non-positive spectral function cannot 
be understood in terms of physical particles on the mass 
shell. On the other hand, a non-trivial quark propagator 
must depend non-trivially on momentum if it is defined 
in the whole complex plane. As will become clear in Sect. 



IXIIII a VMD-based model produces a complex spectral 
function p(u>) on a complex contour C, which results in a 
quark propagator with cuts only. 

In a model without confinement the spectral represen- 
tation of a propagator is a well defined concept. In a 
gauge theory, like QED, the spectral representation de- 
pends on the particular gauge, because the two-point 
function does. In a theory with confinement not much is 



4 



known about the analytic properties of the quark propa- 
gator, except for the fact that poles at real positive val- 
ues of p 2 with positive residues are certainly excluded. 
In QCD such a representation certainly exists in pertur- 
bation theory where confinement is not manifest. The 
study of p(ui) within QCD yields at LO, for m — ► 0, the 
following expression |2Gj: 



p(u>) = 5(lo — m) + sign(o>) 



aC F 1 - £ 
An to 



0(co 2 



,(2.5) 



where £ is the gauge parameter. However, if the gen- 
eral representation is valid, its detailed properties may 
be quite different in the non-perturbative regime. 

As will become clear below, the real strength of the 
ansatz (|l.lf) relies in the fact that having assumed certain 
properties of p(w) reduces the calculations of physical 
observables to nothing more than the standard one-loop 
analysis. In addition, it allows for going from the Eu- 
clidean to the Minkowski space, back and forth. Thus, 
our assumption is essentially that of analyticity of the 
quark propagator S(p), and the possibility to analyti- 
cally continue it in the whole complex plane. This looks, 
in principle, very different from the approach invoked 
in non-local models, where only the Euclidean region is 
used to justify the propagator. Nevertheless, calculation 
of physical observables require in practice an extrapola- 
tion into the complex p 2 -plane, which by itself can only 
be justified through analytic continuation. The previous 
argument does not justify Eq. Ijl.ip. but it shows that we 
are not making any additional assumptions as compared 
to those implied in non-local models. 

Although in perturbation theory the integration con- 
tour may be kept on the real axis, there are cases where 
singularities may pinch the integration path. This cir- 
cumstance becomes a problem, since either analyticity 
or relativistic invariance may be spoiled. In this regard 
a number of prescriptions have been devised in order to 
avoid such a situation |2lL I22I |23| . Thus, in general, we 
will assume that the integration path is an arbitrary con- 
tour chosen in a convenient way. This contour integration 
requires effectively considering complex masses. Genuine 
non-local models formulated in the Minkowski space also 
require a specification of the integration contour in mo- 
mentum space in order to keep relativistic invariance |24| . 
An example of a definite prescription of the choice of such 
a contour is given in Sect. IXIII 

The quark propagator may be parameterized in 
the standard form, 



S{p)=A{p)i> + B(p) = Z{p) 



with 



A(p) = J dio 



j> + Mjp) 

p 2 -M 2 (p)' 
p(co) 



(2.6) 



B(p) 



(III! 



p 2 — Ld 2 ' 

p(lo)lo 

p z — LO z 



(2.7) 



and the mass and wave function renormalization func- 
tions given by 



M(p) 



B(p) 



A{ P y 

Z(p) = (p 2 - M 2 (p))A(p), 



(2.8) 
(2.9) 



respectively. Let us note that if we had p(lo) = p(—u>), 
then the quark mass would vanish, M(p 2 ) = 0, and spon- 
taneous breaking of the chiral symmetry would then be 
precluded. Thus, in general, we expect p{to) not to be an 
even function. In the following sections we will compute 
one-by-one the physical observables and accumulate con- 
ditions that have to be satisfied by the moments of the 
spectral function p(ui). 



B. Quark condensate 

With help of the representation (|l.lfl the quark con- 
densate (for a single flavor) may be straightforwardly 
computed, yielding 



(qq) = -iN c 



-MN r 



dujp(u> 
diop{u>) 



d 4 p 
d 4 p 



Tr 



1 



•p — LO 

LO 



(2tt) 4 p 2 - LO 2 



(2.10) 



where the trace is over the Dirac space, and N c = 3 is 
the number of colors. The integral over the momentum 
p is quadratically divergent. This requires the use of an 
auxiliary regularization method, removed at the end of 
the calculation. With three-dimensional cut-off, A, one 
gets for large A 



W) 



= -4^ JdLOLop(Lo) 



2A 2 + Lo 2 log 



\ 4A 2 



The finiteness of the result at A 
ditions 



Pi = 0, p3 = 0, 



(2.11) 

00 requires the con- 
(2.12) 



and thus 

(qq) - 



N f 

J dcjl °S(^ 2 p 3 P^) 



-^P's- (2-13) 



Exactly the same conclusions are reached if the four- 
momentum auxiliary regularization is introduced. Note 
that the pz = spectral condition allowed for rewriting 
log(a-> 2 /A 2 ) as log(o-' 2 ), hence no scale dependence (no "di- 
mensional transmutation" ) is present in the final expres- 
sion. The dimensional regularization in 4 + e dimensions, 
gives 



dL0p(L0)L0 3 



\0g(L0 2 /p 2 )-- + l 

£ 
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However, the dimensional regularization "hides" some 
conditions, for instance here it leads only to the p 3 = 
condition, and does not require the p\ = condition. 
This is due to the fact that in the dimensional regular- 
ization the power divergences have a fixed ratio. 

Finally, we remark that in the perturbative phase 
with no spontaneous symmetry breaking, where p(to) = 
p(— uj) = S(lo), we have (qq) = 0. With the accepted 
value of 



(qq) = (uu) = (dd) ~ -(243 MeV) 3 



(2.14) 



(at the typical hadronic scale of 0.5 — 1 GeV) [25J we can 
infer the value of the third log-moment, p' 3 . The sign of 
the quark condensate shows that 



Pz >0. 



C. Non-local quark condensate 



(2.15) 



In various QCD studies (see e.g. Refs. |26l l27l l28l l29| 

and references therein) the non-local quark condensate, 



_ (q(0)q(x)) 



(?(0)«(0)>' 
plays an important role. In our approach 

(q(0)q(x)) = -4iN c f dujp(uj) 1 '''' 



(2.16) 



(2tt) 4 p 2 - uj 2 



(2.17) 



and, consequently, after performing the Fourier-Bessel 
transform, 



Q[x) = _ / dup^A^ 1 ^^- 



1 

7 3 



(2.18) 



where x denotes the Minkowski coordinate. A related 
quantity is the average vacuum virtuality of the quarks, 
A g , defined through 



A 



a _ (V(d 2 )<l) 



W) 



(2.19) 



With our spectral regularization method we find that 
(q(d 2 ) n q) fdujp{uj)fd 4 p(-p 2 ) n uj/(p 2 -uj 2 ) 



in particular 



Jdujp(uj)Jd i puj/(p 2 -uj 2 ) 

/ \n Pgn+3 
P3 



(2.20) 



A 2 = (2.21) 

P3 

while the QCD sum rules estimates suggest the value 
A 2 ~ 0.5 ± 0.1 GeV 2 The positivity of A 2 and p' 3 
(see Sect. IIIB|) enforces 



p' 5 < 0. 



(2.22) 



D. Vacuum energy density 

Continuing the quest for the conditions on the spectral 
function p(to) we now study the vacuum energy density. 
The energy-momentum tensor for a purely quark model 
is defined as 



e^(x) 



q(x)- {7^ + Yd 1 "} q(x) - g^C(x). 



(2.23) 



At the one-quark-loop level 



') = -iN c N f / dujp(uj) 



Tr- 



1 



1 



d 4 P 



(7V + 7y)- s '"'()i-w) 



= -AiN c N f 
= Btf" + (6 



du>p(u>) 



d 4 p p^p v -g^(p 2 -lu 2 ) 
(2tt) 4 p 2 -uj 2 

(2.24) 



where Nf is the number of flavors and (6^)0 is the 
energy-momentum tensor for the free theory, i.e. eval- 
uated with p(uj) = S(uj). The quantity B is the vacuum 
energy density given by 



B = -iN c N f J^)J^ ¥2 



(2.25) 



where in the subtraction of the free part we have used the 
spectral condition l|1.2|) . The integral over p is quadrati- 
cally divergent, but there is an additional power of uj as 
compared to the case of the quark condensate. Hence, 
the conditions that have to be fulfilled for B to be finite 
are pi = and P4 = 0. Then 



B 



16tt 



3Nc 
16tt 2 



-Pi 



(2.26) 



for three flavors, Nf = 3, used from now on. 

Interestingly, the even conditions (here quadratic and 
quartic) imply that p(oj) cannot be positive definite] oth- 
erwise the even moments could not vanish. 

An alternative expression for B can be obtained from 
integrating by parts in the variable p and using the con- 
ditions P2 = Pi = 0, 



B = 2iN c N f / dujp(aj 



d A p 



log(p 2 



■ UJ 



(2.27) 



which reminds us of the vacuum energy density calcu- 
lated in the effective action formalism. Another inter- 
esting version of B is obtained by integrating Eq. (|2.25(l 
with respect to the variable po first, when we get 



B 



2N c N f / duop{uj) 



d 3 k 
(2^)3 



Vk 2 



(2.28) 



The interpretation of this equation is obvious: B is a 
weighted integral of negative energy quarks with con- 
stituent mass uj. Upon imposing the pi — and p± — 
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conditions the integration of the three momentum inte- 
gral yields Eq. QTZty . 

In the case of spontaneous chiral symmetry breaking 
one expects, (#oo) < (^oo)o, or B < 0. According to 
the most recent QCD sum rules analysis for charmonium 
[251 Enl . one has, for three flavors, 



-(224±^MeV) 4 . (2.29) 



Note the large uncertainty in this result. The negative 
sign of B enforces 

p' 4 > 0. (2.30) 



E. Low- and high-momentum behavior 



At low momenta we may formally expand the quark 
propagator l|2.2|l to obtain 

S(p) = -±[ d» P M (t) k = -± P+tf&M) 



fc=0 



Hence the low-momentum expansion involves the inverse 
moments l)1.4[l . In particular, for M and Z of Eqs. 
wc fi nd 



M(0) = '— 



dM(p 2 ) 



dp 2 



3 2 =0 



p-l 

P-2 
P-3 
P-2 

P 2 -1 



P-lP-A 

? ' 

P-2 



Z(0) = ^i. 

P-2 



(2.32) 



A knowledge of these quantities, for instance from lat- 
tice calculations, would help to determine the inverse 
moments <|1.4[1 and constrain the spectral function. Ac- 
cording to Eqs. (|2.32l) . the positivity of Z(0) leads to 
p_2 > 0, while combined with the positivity of M(0) 
gives p-i > 0. 

In the limit of large momentum, p — > oo, we formally 
have 



S( P ) 



^ J duip(w) + J dujujp{ 

^ J duLJ 2 p(u) + ... 



(2.33) 



Normalization of the quark propagator in the asymptotic 
region to 1 jj> leads to the condition p — 1. Furthermore, 
since M(p 2 ) should vanish asymptotically, we conclude 
that p\ = 0. 

Note also that if all spectral conditions 211.3(1 are 
assumed then the asymptotic large momentum expan- 
sion would yield a trivial free massless quark propagator. 
Thus, the high-momentum expansion cannot represent 
the full (and non-trivial) propagator. This indicates some 
non-meromorphic structure in p at infinity. In Sects. IXIII 
and IXIIII we will present a particular realization of this 
situation. 



III. GAUGE TECHNIQUE AND THE VERTEX 
FUNCTIONS 



Up to now we have accumulated the lowest spectral 
conditions (|1.2I1.3|) . up to n = 4. Higher conditions will 
appear in the next sections by requesting the twist ex- 
pansion of correlation functions. First, however, we need 
to introduce the coupling of currents to quarks. 

In QCD, the vector and axial currents are defined as 



(3.1) 



J^(x) = q{x)>f^q(x). (3.2) 

Although our formulas below are valid for the SU(Nf) 
flavor symmetry, in this paper we will be concerned 
mainly with pion properties. This will be understood by 
replacing the Gell-Mann matrices, A a , by the Pauli ma- 
trices, T a . Conservation of the vector current (CVC) and 
partial conservation of the axial current (PCAC) implies 
that 



d^\x) = o, 



q(x)M i^B-^l( x )' 



(3.3) 
(3.4) 



with Mo = diag(m„, md, m s ) denoting the quark mass 
matrix. Obviously, any effective theory of QCD must 
incorporate these constraints. CVC and PCAC imply 
a set of flavor-gauge and chiral Ward-Takahashi iden- 
tities among correlation functions involving vector cur- 
rents, axial currents, and quark field operators, which are 
based on the local current-field commutation rules 321, 



J^ a {x),q{x') 



= -l^q{x)S(x - x'), 

x =x' I 

= -75^-q(x)S(x - a/). (3.5) 

Xn=x' n I 



A number of results are then obtained essentially for free. 
In the low-energy regime pions arise as Goldstone bosons, 
and the standard current algebra properties hold. In the 
high-energy regime, parton model features such as scaling 
and the spin- 1/2 nature of hadronic constituents may be 
recovered. If one restricts to the one-quark-loop approx- 
imation, the results also provide a particular solution to 
the large- N c counting rules at the leading order. 

To solve the Ward-Takahashi identities we follow the 
gauge technique proposed in Ref. (see also Ref. 0, 
l33| ) which has the nice feature of linearizing the equa- 
tions, since they deal with unamputated Green functions. 
This is in contrast to the more standard approach of writ- 
ing the Ward-Takahashi identities for amputated Green 
functions, in which case non-linear equations arise. The 
gauge technique has been mostly used in the past as a 
way to obtain solutions to the Schwinger-Dyson equa- 
tions, both in QED [13, H| and in QCD HJ. Only 
recently it has been used to study hadron phenomenol- 

ogy m 
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A. Vertices with One Current 

The vector and axial unamputated vertex functions are 
defined as 

A£V,P) = iS(p')r^ a (p',p)iS(p) (3.6) 
= f d 4 xd 4 x'(0\T{j£ a (0)q(x')q(x)}\Q)e ip '- x '~ ip - x , 
A% a (p',p) = iS(p')Y^ a {p\p)^{p) (3.7) 



d 4 xd 4 x' (0\T {J% a (0)q(x')q(x)} \0)e lp x ~ lp - x , 



respectively. Here the T's represent the corresponding 
amputated vertex functions [67j . The Ward-Takahashi 
identity for the full vector-quark-quark vertex reads 

(p' -p)X'"(p'p) = s(p')y - y s{p) ■ (3 - 8) 

Likewise, for the axial-quark-quark vertex we have 

(p'-p)Mp',p) = 5(p')y75 +75 yS(p) . (3.9) 

The gauge technique, introduced in Ref. [ljj, consists 
of writing a solution for the vector unamputated vertex 
in the form 



<g" -u 1 2 i>-uj' 



(3.10) 



The axial-vertex ansatz reads 



A») = \ dwp^j-^^ 



2uq p ' \ \ a i 



75' 



2 p-uu 
(3.11) 



In this way the Ward-Takahashi identities are linearized. 
It can be readily verified that these ansatze fulfill the 
identities (|3.8I3.9|) . respectively, up to undetermined 
transverse pieces. A consequence of the axial Ward- 
Takahashi identity is the occurrence of a massless pseu- 
doscalar pole identified with the pion, which takes place 
only if p(uj) ^ 6(u>). 

The pion wave function, corresponding to the tt — ► qq 
vertex, is defined as 



K(P + 9,P) = i / d 4 xe- lpx (0\T{q(0)q(x)}\TT a (q)). 



(3.12) 



Near the pion pole we get 



AT(p + q,p) 



M I/'" '!■!'< 2 A °(P + 3'P)' ( 3 - 13 ) 



where the pion wave function is given by 

in/ \ f , , , i 2uj \ a i 

Al(p + q,p)= / dwp(w)-— — -—75- 



ui f„ 2 p 1 - u 



(3.14) 



We recognize in our formulation the Goldberger-Treiman 
relation for quarks; under the spectral integral over u> 
the pseudoscalar coupling of a pion to the quarks is the 
ratio of the spectral quark mass to the pion weak decay 
constant, 



g„(oj) = (3.15) 
3-k 



B. Vertices with two Currents 

The vertices with two currents, axial or vector, will be 
needed below when computing form factors. We define 
the axial-axial vertex (other vertices can be done in a 
similar fashion) by 

(27r)V 4 V + q'- P - q)Ar/\p\ q';p, q) 

d 4 xd 4 x'd 4 y'd 4 ye l{q '- x ' +p '- y '- q - x - p - v) 
x(0\T{j^(x)r/(x')q(y)q(y')} |0), (3.16) 

which in SU(2) fulfills the Ward-Takahashi identity 

-tq»A^ b (p',q';P,q) = i€ bac A^{p' ' ,p) (3.17) 
+y75A^V-9,f) + A^V,p + q)^ 75 . 

Up to transverse pieces one gets the solution 



Aft"; 

AA 



v 2ujq' v \ T b i ( _ u 2ajq^\ r„ 



1 q' 2 ) Tt> '2i> + i-u 

fV II 

ujq' q p \ i 
crossed -I _ n - o a b f- 



2 / ^5^7 
q I 2 



q' z q 2 J p — oj 



(3.18) 



The 7r<7 — > irq unamputated scattering amplitude is de- 
fined as 

(2n) 4 S^(p / + q'-p- q)At(p', q';P, q) (3.19) 
d 4 xd 4 x'e l ( p '- x '~ p - x \lT b {q')\T{q{x)q{x l )}\Tt a {q)). 

At the pion poles, q 2 , q' 2 — > 0, we get 



A^\p',q';p,q) - ^f^{p',q';p,q), (3.20) 

q 2 q 



where 



A 



il(p',q';p,q) = j <^p{u>)-^{^8 ab 

ui i lu •) i . . 

t75T6t _ ~2 "T75 r a + crossed > -j . (3.21) 

/ p + tf-uj / J p- U) 
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IV. VACUUM POLARIZATION 

The vacuum polarization is obtained from the vector- 
vector correlation function, which is constructed by clos- 
ing the quark line in the unamputated vector vertex 
lpTTD|l . with the result 



/IT 

-JV, 



vv (q)= / d±xe-^{0\T{j» a {x)r v b {0)}\Q) 



d 4 p 



-N, 
Tr 



(2^)4 
dujp{uj 



Tr 



d 4 p 



(2n) 4 
A„ i 



:7m- 



2 i>-uP v 2 



_i>- i-u 

We use the dimensional regularization 



n 



vv 



(q) = S ab -g» v + 



q^q v 



(4.1) 
|fi?H and obtain 
U v (q% (4.2) 



with 



IV (g 2 ) = 



N c 



p{uj)du) x 



(4.3) 



-2u> 2 I(q 2 ,Lj) + q 2 (~-I(q 2 ,Cj)) 



where the one-loop functions I(q ,u>) and I{q 2 ,w) are 
introduced in App. [X] We note that the vector wave 
function renormalization, 



n'(o) 



El 

3 



p(u>)du> 



1(0, w) 



(4.4) 



diverges, which is the case of perturbative theories, like 
QED, as well. 

The one loop integral satisfies the twice-subtracted dis- 
persion relation (see App.^J, 



dt lmU v (t) 



t 2 t 



(4.5) 



This is in contrast to non-local quark models formulated 
in the Euclidean space, where the dispersion relation is 
postulated, but never deduced. As a matter of fact, even 
in local models, such as in those with the proper-time 
regularization, dispersion relations do not hold |35| due 
to the presence of essential singularities generating non- 
analytic structure in the complex g 2 -plane. 

To end this section, we compute the cross section for 
the reaction e + e~ — * hadrons. This quantity is propor- 
tional to the imaginary part of the vacuum charge polar- 
ization operator. Asymptotically, at large s, we find 



cr(e + e — > hadrons) 



47r «QED 




where ej is the electric charge of the quark of species 
i. Thus, the proper QCD asymptotic result is obtained 
when the spectral normalization condition (|1.2I) is im- 
posed. 



V. PION WEAK DECAY 

The pion weak-decay constant, defined as 

(0 \JT{x)\ **(?)) - iterate*-*, (5.1) 

can be computed from the axial-axial correlation func- 
tion. We insert a complete set of eigenstates into the 
correlator. 



d i xe~ zqx (0\T{J A ia (x)J A *(0)} |0) 



q^q v 



(5.2) 



and recover the pion pole, with the dots indicating pieces 
regular in the limit q 2 — > 0. The procedure of closing the 
quark line in the unamputated axial vertex (|3.11(l results 
in 



-in 



AA 



(?) 



-JVo 



N c / duip(uj) 



Tr 



ft — <j[ — U 



T 



d 4 k 

W) 
d 4 k 

(2^ 
2ujg f _ l 
q 2 



rTr 



A*X(k + q, k)j„j 5 



h 
2 

(5.3) 



75 



_ „ ^6 

2 ft- w 7 75 2 



Note that the above expression involves only one full 
vertex, (7^ — 2toq fM /q 2 )jc l X a /2, and one bare vertex, 
"fu"f5^b/2. This is needed to avoid double counting, and 
complies to the method of Pagels and Stokar 36]: in a 
two-point correlator all diagrams of the underlying the- 
ory (QCD) can be grouped in such a way as to dress the 
quark propagators (self-energy renormalization), and one 
of the vertices, while the other vertex remains in the form 
from the underlying theory. If both vertices were dressed 
and no additional subtractions were introduced, double 
counting of the diagrams of the underlying theory would 
result. 

With the dimensional regularization Eq. 15.3fl becomes 



Li AA 



(q) = Sab 



-9 



//.;/ 



fl A (q 2 ), (5.4) 



with 



U A (q 2 ) =IL v (q 2 ) + 4N c J duj^pi^Iiq 2 ,^). (5.5) 

As we can see, spontaneous breaking of chiral symme- 
try implies a pole in the axial-axial correlator, with the 
residue proportional to the squared pion weak decay con- 
stant. The result is 

fl = 4A c y dwp(cj)w 2 I(0,u). (5.6) 
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A finite value for /„• requires the condition p 2 — 0. Then 

£ = dulogiu^p^) = -^ P > 2 . (5.7) 

Again, the spectral condition p 2 — guarantees the ab- 
sence of dimensional transmutation for this log moment. 
The value of the pion decay constant can thus be used to 
determine p' 2 . The sign is, obviously, 



P2 < 0. 



(5.8) 



VI. WEINBERG SUM RULES 

The basic idea behind the Weinberg sum rules [37], |3?| 
is that at high energies chiral symmetry breaking should 
be small. There are two equivalent ways to derive expres- 
sions for these sum rules in our model: from the absorp- 
tive parts, or from the dispersive parts of the correlators. 
Both are equivalent due to the dispersion relations, which 
we have shown to hold in Sect. II Vl Here we present the 
derivation from the absorptive parts. The vector- vector 
and axial-axial correlation functions can be subtracted 
from each other, yielding for the imaginary parts 

1 



-{lmU v (q)-lmIL A (q)} 



^ J dwcj 2 p(u)Jl - ^-9(q 2 - 4c 2 ). 



(6.1) 



Next, we integrate with respect to q 2 , and use the spec- 
tral conditions and the definition of / w to get 



1 

N 



da 2 

-\{lmn v (q)-lmn A (q)} 
o 1 



(6.2) 



lim 



2tt 2 a^c 



duiuj 2 p(uj) < lo 



A 2 



f 2 



which coincides with the first Weinberg sum rule. Now, 
if we compute the left-hand side of the second Weinberg 
sum rule, we get 



1 



dq 2 {Imll v (q) -Imfl A (q)} 



AT f Q 

= J£ J dww 4 p(u;)]oga; a = (6.3) 

The result, according to the second Weinberg sum 
rule, should involve on the right-hand side the quantity 
m(qq) — f 2 m 2 , which vanishes in the chiral limit. In- 
stead, our formula involves the vacuum energy density, 
B, which does not vanish. This violation of the second 
Weinberg sum rule is similar to findings in other chiral 
quark models, and reflects, in this regard, a deficiency 
of those models as well as of the present approach (see, 
e.g., the discussion in Ref. Hi|). A study, to be pre- 
sented elsewhere, reveals that this should not be consid- 
ered a drawback of the spectral representation method, 
but rather a feature of the particular solution of the axial 
Ward-Takahashi identity 



VII. PION ELECTROMAGNETIC FORM 
FACTOR 

A. Form factor 

The electromagnetic form factor for a positively 
charged pion, tt + = ud, is defined as 

(7r+(p')|^ m (0)k+(p)) = er c ;») (7.1) 

with q = p' — p. Following the method of Sect. II I II we 
compute the form factor by using the 7rq — > nq scattering 
amplitude, closing the fermion line, and tracing with an 
electromagnetic vertex: 



C(p',p) = -N c 



(2nj 



rTr 



(7.2) 



Through the use of Ea. 13.211 we get 



iW.p) = -iv c / du P (u>) 



Tr 



V2u 
i 



d 4 fc 
(2tt)4 



"75 



75" 



ui ^ + ^ — lu $ — ft — uj 



(7.3) 



For on-shell massless pions the electromagnetic form fac- 
tor reads 



4N, 
I 2 



2 , dwp(tu)uj 2 I(q 2 



(7.4) 



which due to Eq. I|5.6|l is obviously normalized to unity 
at q 2 = 0, F% m (0) = I. With help of App. E|we derive 
the low-momentum expansion, 



F^tf) = 1 



4vr 2 / 2 



q 2 po , q 4 p-2 , q 6 p-4 



60 



240 



(7.5) 



The mean square radius reads 



= 6- 



,dF 
'dq^ 



f=o 4tt 2 / 2 



dcop(uj) 



47T 2 /, 



2 ' 



(7.6) 



which coincides with the unregularizcd-quark-loop re- 
sult [40| and also shows that in the present framework 
the pion is an extended object. The numerical value is 



(r 2 )r| th = 0.34fm 2 



(r 2 )r| oxp = 0.44fm 2 , (7.7) 



which is a reasonable agreement. One should not expect 
a perfect agreement since from the chiral perturbation 
theory it is well known that pion-loop corrections provide 
a sizeable enhancement for (r 2 )" n . We note that the 
knowledge of the pion electromagnetic form factor allows 
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to determine the even negative moments of the spectral 
function, cf. Eq. \7.b} . This will be used in Sect. 1X111 
to build a vector-meson dominance model. Based on the 
properties of the one loop integral I(q 2 ,w) (see App. lA"|l . 
the pion form factor also fulfills a dispersion relation in 
our formalism 



Fr(^) = i + -/ - 

TV In t t 



dt ImF% m (t) 



iW 



(7.8) 



B. Twist expansion and spectral conditions 

In the limit of large momentum we find, according to 
Eq. (S3) , 



Nr. 



,2 /, .2 



dujp(uj)uj 2 {2 log(g 2 ) + 



[log(- g > 2 ) + 1] + ^ 



log(- 9 2 /^ 2 ) - 



(7.9) 



With help of the spectral conditions l|1.3|l for n 
2,4,6, ... we can rewrite this expansion in the form 



K m (q 2 ) 



47T 2 / 2 

8B 



q 2 



g 6 
Q 4 



Q 6 



(7.10) 



with Q 2 = — q 2 . Note the very interesting feature: the 
imposition of the spectral conditions l|1.3l) removed all the 
logarithms of q 2 from the expansion 1(7. 9|1 , leaving a pure 
expansion in inverse powers of q 2 . Thus, factorization 
has been achieved, which in our opinion is one of the 
major successes of the present approach. Conversely, in 
order to obtain factorization, the conditions i|1.3f) must be 
assumed. In the present calculation only even moments 
of the spectral function p(lu) appeared. To involve the 
odd moments one needs to consider a different quantity, 
for instance the scalar pion form factor. We recall that 
the odd spectral conditions were also needed in Sect. Ill Cl 

The leading-twist coefficient in expansion (|7.10|l has a 
very simple physical interpretation: it involves the ra- 
tio of the vacuum energy density, B, and f 2 . Finally, 
we remark that the pure power behavior of Eq. H7.10|) is 
characteristic of a bound-state object, and was obtained 
in non-local models 36] and more recently in the instan- 
ton models |4l| and the Nambu-Jona-Lasinio model Q ■ 

It is worth stressing that all spectral conditions p. 3(1 
arc needed. If we just impose a finite number of them, 
say up to order AT, then there appear logarithmic correc- 
tions starting at order iV+2, of the form \og{Q 2 )/Q 2N+2 . 
This is what happens in the Nambu-Jona-Lasinio model 
when a Pauli-Villars regularization with quadratic sub- 
tractions is used; the pion form factor has proper leading 



twist behavior but a logarithmic contribution at sublead- 
ing twist @. 

Plugging the numbers for B and f n we get for the 
leading twist contribution 



Q 2 K m (Q 2 ) 



= ~— =■ = (0.78 ±0.61) GeV 2 , 

twist — 2 O 

(7.11) 



where the uncertainty in the model value comes from the 
uncertainty in B. The experimental result for the full 
form factor is Q 2 F^Q 2 ) = 0.38 ± 0.04GeV 2 as taken 
averaging some old ^2 ar± d recent data as compiled 
in Ref. in the region 2GeV 2 < Q 2 < 6GeV 2 . On the 
other hand, a remarkable finding from the data is that 
a vector meson dominance monopole model obeying the 
pion charge radius up to Q 2 ~ 1.6GeV 2 |45j . 



FTiQ 2 ) 



A 2 



Q 2 + A 2 



(7.12) 



This yields Q 2 F° m (Q 2 ) = 0.41 - 0.45 GeV 2 depending 
on whether one takes A = (6/(r 2 )) 1 / 2 = 0.73 GeV or 
A = My = 0.77 GeV, respectively. This corresponds to 

2 F™ 1 (Q 2 ) = A 2 = 0.53-0.59 GeV 2 for the same 

twist — 2 

range of A values, to be compared with our estimate, 
Eq. ifTTTJ) . Motivated by the VMD success in describing 
the pion form factor in the space-like region, we will study 
further consequences of this scheme in Sec. IXIII 



VIII. ANOMALOUS FORM FACTOR 

A. Vertex function 

The AVV- vertex is defined as 
(2tt) 4 (S(p -qi + q 2 )A^y' a ' ' b (p, q 1 ,q 2 ) 
= i J d i x\d A X2d i x x (8.1) 

{0\T{j^ c (x)J^ a ( Xl )J^\x 2 )]\0)e l ^ x -^ + ^ X2 \ 

where q\ is ingoing, while p and q 2 are outgoing. The 
solution fulfilling the relevant Ward-Takahashi identities 
can be written. Going to the pion pole, p 2 — > 0, yields 

b P^ 1 b 

A^y' a ' ,(3 - h (p, qi,q 2 ) -> -^T^ ;/3 ' 6 (p, qi,q 2 ). (8.2) 

For a neutral pion, 7r°, and two photons one gets 



O (9i,g2) = -N c / dujp(w) 



Tr 



CO 

■~ri5To- 



crossed, 



d A k 
iQY- 



(8.3) 
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where Q = B/2 + I 3 = 1/2N C + r 3 /2 is the quark charge 
operator. Straightforward calculation of the traces yields 

r^o 77 (p,<7i,<72) = e^apq^p 13 Fn 1 ~ l (p,q 1 ,q 2 ), (8.4) 

where the pion transition form factor, 

8 f 

F-K-rt{p,Qi,<l2) =-j- J dwp(u)uj 2 K(p 2 ,ql,ql,u), (8.5) 

has been introduced, and the three-point loop function, 
AT, is presented in App. [S] 

B. Neutral pion decay 

The amplitude for the neutral pion decay, 7r°(p) — » 
7(^1, p) + 7(92, v), can be directly computed from the 
former expression by taking the on-shcll photons, q\ — 
q\ = 0, and the soft pion condition, p 2 — 0. We find 

1 f d 4 k 1 



2^(0, 0,0) = -- ] d.p^u . j (27r)4(fc2 _ w2) 3 



^ 2 U 



du>p(u>) 



^ 2 f, 



(8.6) 



which, when the spectral condition (|1.2|) is used, coin- 
cides with the standard result expected from the QCD 
chiral anomaly. 



C. Transition form factor 

For two off-shell photons with momenta q\ and qi it is 
convenient to define the photon asymmetry, A, and the 
total virtuality, Q 2 , 



A = 
Q 2 ■ 

or, equivalently, 



? 2 

i\ + <i ' 



-1< A < 1 



•7) 



{ ^Aq 2 , 



At the soft pion point we find 



2r(o,- 



dwp(w)w 2 x 



(8.9) 

q 2 ,^). 



Through the use of expansion l|A18|l and the spectral 
conditions l|1.3[l for n = 2,4,6,... we may write, after 
straightforward manipulations, 



^*7T7*7*(Q 5^) — 

2p 2 



1 



da; x 



^(l + ^ 2 (2x-l) 2 ) 



(8.10) 



Q 2 (l-A 2 (2x-1) 2 ) Q 4 (l - A 2 (2x - l) 2 ) 



We can now confront this expression to the standard 
twist decomposition of the pion transition form factor 

m, 



which via Eqs. (|5.7l2.26jl yields 
4/. '•' 



J^(A) 
JW(A) 



N c jo 
8/^A 2 



d.r 



V?i 2) (x) 



1 - (2x - 1)M 2 ' 

1 , ^i 4) (x)[l + (2x-l) 2 v4 2 i 
dx 



.11) 



.12) 







[1 - (2x - 1) 2 A 2 Y 



with 



8B 



(8.13) 



(8.14) 



Note that this is exactly the same combination of B and 
fn as in Eq. I|7.11|l . Numerically, we get 



A 2 = (0.78 ±0.61) GeV 2 



(8.15) 



An estimate made in a non-local quark model of Ref. [42i 
E| provides A 2 = 0.29 GeV 2 . 

The form of the expansion coefficients (|8.12I8.13|) 
shows that the twist-2 and twist-4 distribution ampli- 
tudes for the pion are, at the model working scale Qo, 
constant and equal to unity. 

It is interesting to look at higher-order twist coeffi- 
cients. With help of App. ^ and with the spectral con- 
ditions (|1.3fl we get 



n+2 



X 



1 



1 



/2-n) 

n+1 

, , ■ (8.16) 



1 + A(2u-1)Q 
This yields the result 

tpl n) {x) = 6(x)9(l-x) for n = 2,4,6, 



.17) 



All these amplitudes are by convention normalized to 
unity. The prediction of the model is that they do not 
depend on the Bjorken x variable at the scale Qq. 
In the limit q\ — q%, or A = 0, we get 



-^7*7*7r(Q : 0) 



4/* 

n c q 2 



1 



2A 2 



(8.18) 



An analogous calculation to the one presented in 
Refs. |8|, lUJ produces the following light-cone pion wave 
function in the present model: 



#(x,/c_l) 



4vr3/ 2 



du>p{oj) 



-6{x)6{l-x). 



(8.19) 



12 



Again, this form corresponds to the low-energy scale of 
the model. The pion light-cone wave function H8.19fl sat- 
isfies at k± = the following condition, 

0) = ^iWO, 0, 0) = (8.20) 

In QCD one has a similar relation holding for quantities 
integrated over x I n our model this is inessential 

due to the fact that the x-dependence is constant. This 
triple identity, although comes out easily in our model, 
is difficult to get in local chiral models since there is a 
conflict between recovering the proper anomaly and ob- 
taining factorization of the form factor at high photon 
virtualities (see a detailed discussion in Ref. 

D. QCD evolution 



we may compute the distribution amplitude for any 
value of Q 2 . The result extracted in Ref. [53 and con- 
firmed in Ref. .51] from experimental data |52j| provides 
02(2. 4GeV) = 0.12 ±0.03, hence we can fix the evolution 
ratio to the value 

a(Q = 2.4GeV)/a(Q ) = 0.15 ± 0.06, (8.26) 

which reproduces 02 obtained in our model. Then we can 
predict 

a 4 (2.4GeV) = 0.06 ±0.02 (exp :- 0.14 ± 0.03 T 0.09), 
a 6 (2.4GeV) = 0.02 ±0.01. (8.27) 

The overall picture at the leading-twist and with leading- 
order QCD evolution is very encouraging. For further 
details the reader is referred to Ref. Il7l. 



The results of the previous section referred to the soft 
energy scale of the model. In order to compare to exper- 
imental results, obtained at large scales, the QCD evo- 
lution must be performed. The procedure has been dis- 
cussed in detail in Ref. ^3] , hence here we only sketch the 
method and mention, for completeness, the most impor- 
tant outcomes. For the twist-2 pion distribution ampli- 
tude the leading-order QCD evolution is made in terms of 
the Gegenbauer polynomials. One begins by interpreting 
our result as the initial condition, 



<^ (2) (x,Q ) =6{x)6{l-x). 
Then the evolved distribution amplitude reads 



.21) 
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(2) (*,Q 2 ) 



&x{l^x)Y J Cl ,2 {2x-l)a n {Q), 



n=0 



.22) 



with 



a n (Q) 



2n + 3 



a(Q 2 ) 

3(n + l)(n + 2) \a{Q 2 ) 



7i°V(2/3o) 



x / dxC z J 2 {2x- (x,Q 2 ), (8.23) 



where C 



.3/2 



are the Gegenbauer polynomials, and 

n n+l 



In 3 



(n + l)(n- 



2 ) k k 



Pa = yJV c -|jV>=9. 

With our initial amplitude we immediately get 



(8.24) 



IX. THE 7 -> 7r+7r°7r- DECAY 

In this section we consider an example of a low-energy 
process involving a quark box diagram. The amplitude 
for the decay of the photon of momentum q and po- 
larization e into three pions of momenta pi, 7(9, e) — > 
7r + (pi)7r°(p 2 )7r~(P3), is equal to 



T- 



7(9, e)— ►7r+(pi)7r (p 2 )-T~ (P3) 



= 6i 



dujp{uj) 



Tr 



2NJ ^-i>i-h-h-w V .U 



W 



fn J f - i>\ - W \ frr 



U) 

T~T 



ft — W 



F(P1,P2, P3)£ a f3*re a PiP2Pl . 



(9.1) 



In the limit of all momenta going to zero we get, with 
the condition i|1.2|) . 



F(0,0,0) = 



1 



dajp(u>) 



1 



(9.2) 



dxCl ,2 {2x~ l)tpV\x,Ql) = 1. (8.25) ( 



What actually matters in this analysis is the evolution 
ratio a(Q 2 )/a(Q 2 ). With help of Eq. i8. 2^18.2^18.2^1) 



4tt 2 /J J —rx-j 47r2/ s> 
which is the correct result [53l IH4I l55l | . 

X. PION STRUCTURE FUNCTION 



As we have said in the Introduction, one of the ad- 
vantages of our model over other formulations is that 
calculations can be undertaken both in Minkowski and 
Euclidean space. This proves crucial in the calculation 
of the pion structure function. We recall here that a 
Euclidean formulation allows only for the calculation of 
a finite number of moments of the structure function, 
x n ), for integer n, requiring a subsequent reconstruc- 
tion of the distribution function. This is a cumbersome 
situation. Experimental data are directly obtained in the 
x-space, but structure functions are difficult to pin down 
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for large values of x (typically x > 0.65). That means 
systematic uncertainties for higher-order moments. 

Our following calculation also illustrates an interest- 
ing point. In the Bjorken limit it is assumed that inte- 
grals are convergent fast enough to allow to convert the 
forward Compton amplitude to a quark-target scattering 
amplitude |l4| . We note here that while the former corre- 
sponds to a closed quark line, the latter refers to a quark 
propagator, i.e. an open quark line. In local models, 
such as the Nambu-Jona-Lasinio model, the difference 
becomes subtle (see, e.g., the discussion in Ref. @|) be- 
cause it is not obvious how to regulate open quark lines. 
As we show below, a rewarding feature of the present ap- 
proach is that the connection from the forward Compton 
scattering amplitude to the quark-target scattering for- 
mula prevails, due to the the spectral regularization of 
the vertex functions. As a result, the relation between 
gauge invariance and proper normalization of the PDF 
remains valid. 



A. Derivation from the forward Compton 
amplitude 

The hadronic tensor for inclusive electroproduction on 
the pion reads 



W^(p,q) = —lmT^ v (p,q) 



= W 1 (q 2 ,p-q) [-9^ + 



(10.1) 



W 2 {q\p-q) 



p ■ q 



Pv 



p-q 
—o-q v 



where the forward virtual Compton scattering amplitude 
on the pion is defined as 

T ltu (p,q)=i J d 4 xe^^7r(p)\T{j™(x)jr(0)}\7r(p)). 

(10.2) 

The amplitude can be obtained by taking the residue of 
the double pion pole in the AV — > AV amplitude. The 
gauge invariance requires considering not only the box- 
like diagrams but also the process irj — > tt — > ttj. This 
process may be relevant at low energies, however it does 
not contribute at high energies, since it provides higher 
twist contributions. The hand-bag diagrams yield 



iT„„(p, q) 



Tr 



1 



lST a - 



-N, 



du>p(u>) 



d 4 k 



UJ ft + $ — UJ 



Qlv 



1 



-75T&- 



+crossed. 



— CjJ p 1 — $ — UJ 

(10.3) 



In the Bjorken limit we can make the customary approx- 
imation, 



1 



ft + </i - w 



7 



Lt V 

T 77 



(10.4) 



2 [5^ Q/? 7/5 + ^^7/375j 



c livaf3 r 



The hadronic tensor is obtained as the imaginary part 
in the (p + q) 2 channel. The Cutkosky rules amount to 
making the replacements 

-^—z (-2m)0(p o )S(p 2 - uj 2 ) = (-2m)S + (p,uj), 



hence 



W^ v = -N c J dwp(uj)(j- 



(10.5) 



(10.6) 



d*k 

WT 2 



(-2m) 



,5+(p- k,uj)5 + (q + k,uj) 



(k 2 



,2\2 



Tr 



lb T a (ft + uj)Q^(ft + i + uj)Q^(ft + uj)x 



ihn(j/> - i - u)] . 

The calculation of the traces is straightforward and the 
Bjorken limit of the discontinuity can be found in Adp.IbI 
The result is rather simple, 



W„ v (p, q) 
F(x) 



1 



ImT^(p, q) 



(10.7) 



9iw+ q2 q 2\P» 2x KPv 



2x' 



with 



F ( x ) = l ^(x) + q t (a 



i—u,d.s 



We take tt + for definitcncss and get, 

u n (x) = dn(l — x) = 6(x)6(l — x), 



(10.8) 



(10.9) 



independent of the spectral function p(u>). Thus we have 
recovered scaling in the Bjorken limit, the Callan-Gross 
relation, the proper support, and the correct normaliza- 
tion. This is the same result as found by one of us 
when computing the structure function from the forward 
quark-pion scattering amplitude, Eq. I|3.21|l . in the light- 
cone coordinates. The result has also been obtained pre- 
viously by several means within the Nambu-Jona-Lasinio 
model El Hi 

As one can see, the normalization integral coincides 
with the normalization of the pion electromagnetic form 
factor, 



dxq(x) 



dxq(x) = F° m (0) = 1. (10.10) 
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In addition, we have the crossing property, 

q(x)=q(l-x). (10.11) 

The /c^-unintegrated parton distribution can be shown 
to be equal to 



q(x,k±) = A _^ n I dujp(uj)--^ — -9(x)0(l - x), 



47T3/2 



(10.12) 



which is the same form as in Eq. H8.19fl . hence at the 
working scale of the model, Qo, one has the interesting 
relation 



in the model. As it is well known [Til \El\ . the matrix 
element of the energy momentum tensor coincides with 
the first moment of the PDF 



dx xq(x) = / dxxq(x) = — . (10.19) 
Jo 2 



Actually, this property is a simple consequence of the 
crossing property q(x) = q{l — x) and the normalization 
condition. 



C. QCD evolution 



q(x, kj_) = q(l - x, k±) = $>(x, k±). (10.13) 

valid in our model in the chiral limit. A similar iden- 
tity has also been found in the Nambu-Jona-Lasinio 
model jjjj- Combining Eq. (|10.13|1 with the anomaly 
condition l|8.20[l . we get the following normalization for 
the unintegrated parton distribution at k±_ = 0, 



Q(x, 0_l) 



(10.14) 



Finally, via integrating with respect to k± the following 
identity between the PDF and the PDA is obtained at 
the scale Qo: 



q(x) = ip n (x). 



(10.15) 



This relation holds also in the Nambu-Jona-Lasinio 
model with the Pauli-Villars regularization |l7j . 



The QCD evolution of the constant pion structure 
function has been treated in detail in previous works 
|15l Il6| at LO and NLO order. Nevertheless, in or- 
der to make the paper more self-contained we present 
here the main points arising from that discussion. The 
previous subsections yield the form of the leading-twist 
contribution to the pion structure function at a given 
renormalization point, Qo- The QCD radiative correc- 
tions generate logarithmic scaling violations, which can 
be included in our model by the DGLAP equations [Hsj . 
In particular, the non-singlet contribution to the energy 
momentum tensor evolves as 



Jdxxq(x,Q) ( a(Q) 



Jdxxq(x,Q ) \a(Qo) 



7{ 0) /(2/3o) 



(10.20) 



B. Momentum sum rule 



where 7^ and /3q are given in Eq. H8.24JI . In Ref. it 
was found that at Q 2 — 4GeV 2 the valence quarks carry 
47 ± 0.02% of the total momentum fraction in the pion. 
Downward LO evolution yields that at the scale 



In our formalism the pion expectation value of the 
energy-momentum tensor l|2.23l ) is 



Qo = 313+^MeV 



(n a (qW(0)\*"(q)) - -N c I dup{u) [ — l x 



d 4 k 



rTr 



75 T a 



ft + <j[ — cj 5 b 'ft — 



(10.16) 



i(^ + ^V*)- fl ^- w )|-L- 



Though the use of the spectral conditions we get, for 
on-shell massless pions, 



(n a (q)\en^(q)) = 2 



2 B 
N~fJ! 



The connected piece becomes 

(K a {qW v \-K b {q))c=2q>*q v . 



<5 a6 .(10.17) 



(10.18) 



Thus, the quarks carry all momentum of the pion, as it 
should be, since there are no other degrees of freedom 



(10.21) 



the quarks carry 100% of the momentum. The agreement 
of the evolved PDF 0, 0] with the data analysis [5j| 
is quite impressive. Equation (|10.15f) has been shown in 
Ref. 0| to produce a very interesting integral equation 
relating the evolved PDF and PDA, valid at the leading 
order QCD evolution. 



XI. GASSER-LEUTWYLER COEFFICIENTS 

The gauge technique provides a way to deal with open 
quark lines, an advantage over traditional chiral quark 
models, but it is unnecessarily complicated when dealing 
with processes with closed quark lines. For such a situa- 
tion the effective action approach provides a much more 
efficient calculational tool. It also yields a closer connec- 
tion to previous approaches such as bosonized versions of 
the Nambu-Jona-Lasinio model. The one-quark-loop ef- 
fective action that incorporates the quark-pion coupling 
obeying the Goldberger-Treiman relation l|3.14|l can be 
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written in the form 

S = —iN c J d 4 x J dujp{uj) x 

Trlog [ifi - uj exp (i"fyT a (j) a (x)/ f^)) ■ (11.1) 

This form is manifestly chirally symmetric, with <j> de- 
noting the non-linearly realized pion field. Note the 
formal similarity with a generalized Pauli-Villars regu- 
lator. One may evaluate the Gasser-Leutwyler coeffi- 
cients 0| through the use of standard derivative 
expansion techniques |62| in a similar fashion as done 
in Ref. [UlijJ]. With the lfT~^|) condition imposed, the 
calculation is equivalent to standard quark-model calcu- 
lations with the cut-off removed. The resulting values of 
the Gasser-Leutwyler coefficients are 



h = -N c 
h = N c . 



(11.2) 



Other low energy constants, such as ^3 and I4 require a 
specification of explicit chiral symmetry breaking within 
the quark model. External gauge fields may also be cou- 
pled resulting in predictions for l 5 and l 6 , although any 
choice reflects a particular selection of transverse pieces. 
This and related issues are postponed for future studies. 
It is nevertheless interesting to anticipate here a dimen- 
sional argument which shows why the relevant spectral 
condition for the terms involving fourth-order derivatives 
is, precisely, po = J dujp(uj) = 1. In the case of the 
7T7T scattering in the chiral limit, we have a box diagram 
with four quark propagators, i/(j/) — u>), and four external 
pion lines, each contributing a factor of uj/ f^, due to the 
Goldberger-Treiman relation (|3.15[l . If we are after the 
coefficient with four derivatives we need four additional 
powers of momenta in the denominator, which we may 
account for by squaring the fermion propagator. Thus, 
in obvious dimensional notation we have after adjusting 
the dimensions, 



p(uj)duj 



d A p ( uj 



(2^) 4 V/. 



2 2 

p A — UJ Z 



1 



p(w)dw- I (a^) 4 

J TV 



(11.3) 



due to the fact that the dimensions of the convergent 
integral are set by the spectral mass uj. This shows that 
the terms of dimension four in the effective Lagrangian 
are proportional to pQ. 



XII. VECTOR-MESON DOMINANCE 

Up to now, our considerations have been made for a 
general spectral function fulfilling a set of properties re- 
garding their moments and log-moments. It is quite nat- 
ural to ask whether such a function exists and what are 
the phenomenological consequences of making specific 
ansatze for this function. In this section we construct 



explicitly the spectral function using the phenomenolog- 
ical guidance of the previous sections. 

Some interesting consequences and insight may be ob- 
tained in the present chiral quark model if the vector- 
meson dominance of the pion form factor is assumed, 



Fy(t) 



M 2 
M 2 +t 



(12.1) 



with My denoting the p-meson mass. This form fits the 
recent data |43( remarkably well. As will be shown below, 
the model for the spectral function becomes explicit and 
further interesting results may be obtained. 



A. Vector-meson dominance in the spectral 
approach 

The vector form factor obtained in Eq. (|7.4|) reads, 
through the use of the Feynman parameterization (|A4|) , 



Fy(t) = FT(t) = 



(12.2) 



. . dujp(uj)uj 2 I dxlog [uj 2 + x(l - x)t] . 

47r fir J JO 

If we make a series expansion in t, the integral in x can 
be carried out order by order, hence 



F v (t) 



dujp{ui)uj' 1 



* / dx[x(l-x)] n ± — - 1 



(12.3) 



Nr. ^ 2- 2 "- 1 0Fl> + l) , 

2 , P2-2n „w„ 1 o >n\ 



nT(n + 3/2) 



As one can see, in our model the pion form factor is 
encoded in the negative even moments. Through vec- 
tor meson dominance we get immediately, by comparing 
Eq. (|12.3|) to the expansion of Eq. (|12.1|) . the following 
identification 



P2-2fc 



2 2fc+3^3/ 2/ 2 fcr (fc + 3/2) 

N c M v k T(k + 1) ' 
k = 1,2,3,... (12.4) 



In particular, the normalization condition, p = 1, yields 



Ml = 



Nr 



(12.5) 



This relation is usually obtained when matching chiral 
quark models to the vector-meson dominance and ap- 
pears all over the literature, yielding a quite reasonable 
estimate for the p meson mass, My — 826 MeV with 
/„ = 93 MeV, and My = 764 MeV with f n = 86 MeV 
in the chiral limit. 

The interesting and remarkable point about Eq. (|12.4|) 
is that even though we have determined the negative 
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even moments of the spectral function, the positive even 
moments, obtained by analytic continuation in the in- 
dex n, unexpectedly but most desirably, fulfill the spec- 
tral conditions 1)1. 3[) for the positive moments due to the 
fact that r(n) has single poles at non-positive integers, 
n = 0, —1, —2, . . . Hence 



p2r. 



0, 



1,2,3. 



(12.6) 



Thus, it makes sense to evaluate the log-moments 111.511 . 
since the absence of dimensional transmutation is guar- 
anteed. The log-moments are most easily evaluated by 
analytically continuing the moments to the complex n- 
plane. We then have 



Pn = 



dumj n \og(ijj )p(uj) — 2— / dujuo z p(uj) 



and 



« d 
2—p z 

dz 



P2n = 



(12.7) 



M 2 



r(n)r(|-n) 



r(|) 

= 1,2,3. 



(12. 



where we have used Eq. (|12.5|l . The first few values are 

4/ 2 7T 2 



P 2 



Pa 



N c ' 
2/ 2 Af 2 ^ 2 

N c ' 
2f 2 M v n* 

N r . 



(12.9) 



Since p' 2 and p'^ determine and B, respectively, see 
Eq. I|5. 712.26(1 , we may write the following interesting re- 
lation coming from the vector-dominance model and the 
spectral approach: 



B = — 



N C M V 



N c 64tt 2 
= -(202- 217 MeV) 4 . 



(12.10) 



The uncertainty stems only from using either or My as 
input. Our value agrees within errors with the estimate 



B. Inverse problem 

Although for practical calculations the moments seem 
to contain the relevant information that can be used in 
practical applications, it is nevertheless very interesting 
to write down an explicit formula for the spectral func- 
tion. The mathematical problem is then to invert the 
formula 



P2n 



dujuj 2n pv(u), 



(12.11) 



where the moments are given by Eq. I|12.4(l . The solu- 
tion to the problem is given by the following surprisingly 
simple function, 



< \ 11 



1 



2wioj (1- Auo 2 /M v ) d v 



with 



d v = 5/2, 



at which case we have 



1 3tt'M v ^ 1 1 
2ni AN C w (M 2 /4- w 2 ) 5 / 2 ' 



(12.12) 



(12.13) 



(12.14) 



The function pv(uj) has a single pole at the origin and 
branch cuts starting at ± half the meson mass, lj = 
±My/2. The contour for computing the spectral mo- 
ments is depicted in Fig. ^ The contributions encircling 
the branch points cancel provided dy is half-integer. The 
contribution at infinity cancels, and only the residues at 
the origin contribute. One may explicitly verify with no 
difficulty that the use of Eq. I|12.12|) in Eq. (|12.11f> with 
the contour of Fig. ^ reproduces Eq. I|12.4fl • 

With help of the explicit formula l|12.12(l several inter- 
esting features may be pointed out. The spectral func- 
tion is genuinely complex and it is defined on a complex 
contour C. This precludes positivity conditions |70j . 

It is also interesting to note that in the limit My — > oo 
one gets 



p v {u) 



1 1 

2ni io ' 



(12.15) 



i.e. the massless free theory. This suggests a multiplica- 
tive effect of chiral symmetry breaking on the spectral 
function. Finally, Eq. (|12.12l) cannot be interpreted as a 
constituent model for which one essentially has a pole 
at lu — M in the complex plane (or, equivalently, a 
5(u> — M)), with M denoting the constituent quark mass. 



XIII. QUARK PROPAGATOR IN THE MESON 
DOMINANCE MODEL 

A. Scalar spectral function 

In the construction of the vector spectral function we 
have used the vector-meson dominance principle, which 
has a firm phenomcnological justification. For the case 
of the scalar spectral function we will proceed differently, 
more heuristically. First, we will propose its form in an 
analogy to the form of py. Then, in Sect. IXIII CI we will 
confront our hypothesis to the recent lattice data on the 
quark propagator 0, 0] . 

The scalar function has to satisfy the conditions i|1.3|) 
at odd positive values of n. The analysis of the previous 
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-M/2 



w-Complex Plane 



M/2 



and for the positive even and odd log-moments, 



Mg-V T(d v - n)r(n) 
4 



r(d v ) ' 

M|\ n - 2 47r 2 (OT)r(ds-n)r(n) 



(13.5) 



/>2r, 



P2 ™- 1 ! 4 y 7V C r(d s -2) ' 
n= 1,2,3,... 

The value dy = 5/2 should be used for the vector domi- 
nance model. 



B. Quark propagator 



FIG. 1: The contour in the complex tj-plane for the spectral 
functions in the meson dominance model. The quantity M 
denotes the generic meson mass, My for the vector channel, 
and Ms for the scalar channel. The branch cuts start at 
uj — ±M/2. The pole at ui — occurs for the vector case 
only. 



section suggests the following form 



1 16(d s -l)(d s -2)p' 3 
2m M%(1 - 4w 2 /Mt) d s 



(13.1) 



where the normalization is chosen in such a way that 
the third log moment, p' 3 — —Air 2 (qq)/N c , is reproduced. 
In other words, we fix the normalization with the quark 
condensate. The admissible values of d$ are half-integer, 
since only then the integration around the half-circles at 
the branch points in Fig. ^vanishes. In Sect. IXIII CI the 
preferred value will turn out to be 



d s = 5/2. 



(13.2) 



One may verify that the integration with the prescrip- 
tion of Fig. ^ yields 



P2fe-i = / dcuuj 2k VsM 



(13.3) 



4 3 - k ir 2 (qq) T(d s -k) 
~ ~ Mg 2k+A N C T(d s - 2)r(l-fc)' 

and conditions ((1.311 become satisfied for odd positive n = 
2k — 1. The analytic structure of ps(w) is similar to the 
case of py(ui), except for the absence of the pole at u> = 0. 

For completeness, we list the result for even and odd 
negative moments, 



P-2n 



P-1n-\ — 



M 2 



T(dy + n) 
T(d v )T(n + l) 



(13.4) 



M 2 



47T 2 (gg) r(d a +7i) 
N c F{d s - 2)r(n + 1) 



0,1,2,. 



A straightforward calculation with Eq. (|12. 12113.1(1 

yields the A(p 2 ) and B(p 2 ) functions of Eq. 1(2. 7JL namely 



A{ P 2 ) 
B(p 2 ) 



{1-Ap 2 /M 2 ) d v 
U{d s ~2)(d s ~lW(qq) 



Ap 2 /M 2 ) 



(13.6) 



We note that the apparent pole in A(p 2 ) is canceled when 
the expression in brackets is expanded, and both func- 
tions have no poles in the whole complex plane. The 
functions 1(13.6(1 have branch cuts starting at p 2 = AM 2 , 
where M is the relevant mass. 

The absence of poles, achieved in a rather natural fash- 
ion in our approach, is very appealing, but not completely 
surprising a posteriori. In local chiral quark models, 
where the propagator is usually assumed to be a mero- 
morphic function (with a pole at the constituent quark 
mass), meson vertex functions naturally inherit the dis- 
continuity structure implied by the Cutkosky rules and 
unitarity. In our case, it is the meson form factor which 
is taken to be a meromorphic function through the VMD 
model; unavoidably, the quark propagator must have no 
poles and a certain cut structure conspiring with the uni- 
tarity at the one-loop level in order to produce such a 
form factor with no cuts. 

Alternatively, instead of A and B one may consider 
the more customary mass function, Al(p 2 ), and the 
wave function renormalization, Z(p 2 ), given by Eq. ((2.6(1 . 
They can be written as 



M(p 2 
M 



4d v p 2 



( Ml y 



M 2 



( mi y 



- 1 



Z(p 2 ) = 1 



My 



My ~ A P 2 



(13.7) 



(13.8) 



l&d 2 v M{0) 2 p 2 



( m% \ 



2d s 



My 



( mi V 
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where Mo = M(0) is the value of the mass at the origin. 
We find 

_ l6(d s ~l)(d s ~2)M^(qq) 
Mq — 



d v M^N c 



Z(0) 



AdyM{Qf 
M 2 



(13.9) 



At high Euclidean momenta, Q 2 = —p 2 — > oo, we obtain 
M(Q 2 ) 



i2 _ d v M Ml ( Atp ds 1 



M 2 \4Q 

*M">-»-(^ 



4d 2 v M(0) 2 M 2 ( M|^ 2rf5 " 3 



4Q 2 



(13.10) 



We note that for half-integer ds the tail of M contains 
odd powers of Q, 



M(Q 2 ) 



1 



(<3 2 ) 



2Ws-l 



(13.11) 



and for d s = 5/2 drops as l/Q 3 , for d s — 7/2 as l/Q 5 , 
eic. The wave-function normalization, Z(Q), has the cor- 
rect asymptotic behavior, Z(Q) ~ 1. 

Given the cut structure of the functions A(p 2 ) and 
B(p 2 ) we may look back at the high energy expansion 
(|2.33|) . According to the spectral conditions (|1.3fl one 
would deduce from the high energy behavior that the 
full propagator coincides with the free one. The puz- 
zle is resolved by realizing that the branch cut running 
from ±M/2 to ±oo implies a fractional power behav- 
ior, and hence the function cannot be represented by a 
power series expansion around infinity. From this point 
of view Eq. H2.33f) just expresses the fact that the integer 
power coefficients are exactly zero, as clearly follows from 
Eq. jnSJ. 



C. Comparison to the lattice data 



The decrease of the quark mass as l/Q 3 at large 
Euclidean momenta is favored by recent lattice calcu- 
lations [Tsl llfij . where the fit to the functional form 
c /((Q 2 ) b + (A 2 ) b ) is best when the parameter b is close to 
3 pf. We have performed a x 2 fit for M (Q) of Eq. (|B7|) 
to the data of Ref. 0, . These data, for the case of 
M, have been extrapolated to the chiral limit. We have 
treated Ms and Mq as free parameters. The x 2 fit results 
in the following optimum values, 



M = 
M s = 



303 ± 24 MeV, 
970 ±21 MeV, 



(13.12) 



with the optimum value of \ 2 per degree of freedom equal 
to 0.72. The corresponding value of the quark condensate 
is 




N 
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FIG. 2: The dependence of the quark mass M (top) and 
the wave function normalization Z (bottom) on the Eu- 
clidean momentum Q obtained from the meson dominance 
model with ds = 5/2. The lattice data, extrapolated to 
the chiral limit for the case of M, and taken at the cur- 
rent quark mass ~ 59 MeV for the case of Z, are taken from 
Ref. |isl rTgjl . The thicker lines correspond to optimum pa- 
rameters of Eq. 113. 121 . while the thin lines indicate the un- 
certainty at the one-standard-deviation level. 



-(243.0tg;8 MeV) 3 . 



The functions M(Q) and Z(Q), evaluated at optimum 
parameters, are shown in Fig. [5] with thick lines. The 
thin lines indicate the uncertainty at the one-standard- 
deviation level. The agreement with the data is very good 
for the case of M, and the l/Q 3 fall-off is clearly seen. In 
fact, fitting of the model with values of ds higher than 
5/2, which results in faster asymptotic decrease, results 
in a much worse agreement with the data. As seen from 
the bottom part of Fig. [21 for the case of Z the agree- 
ment is not very good, but we should keep in mind the 
simplicity of the present model and the freedom in the 
scalar channel. For instance, the scalar spectral function 
can be multiplied, without loosing any of the general re- 
quirements, by an entire function. We also wish to stress 
that the optimum value of (qq) obtained by fitting the 
model formulas to the lattice data agrees with the esti- 
mate <|13.13[1 . 

It is interesting to note that even though the mass 
function, M(p ), presents a pole for time-like momenta, 
i.e. there exists a solution to the equation M(p 2 ) —p 2 = 
(13.13) 0, it does not correspond to a physical particle. This 
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Q 2 /M 2 (0) 




-0.1 



-0.08 ^Oq§ ^004 -0.02 



0.02 



Q 2 [GeV 2 ] 



FIG. 3: Square of the quark mass M 2 (solid line), minus 
the square of the Euclidean momentum Q 2 (dotted line), and 
the wave function normalization Z (dashed line), plotted as 
functions of the square of Euclidean momentum, Q 2 . The 
arrow indicates the point where M 2 = —Q 2 , but where also 
Z — 0. As the result, the quark propagator has no pole. From 
the analyticity properties of A and B it follows that the quark 
propagator has no poles in the whole complex-Q 2 plane. 



is because the normalization Z(p 2 ) also vanishes for the 
same value of p 2 . This in fact is just a manifestation 
of the analyticity properties of A(p) and B(p) discussed 
above. Figure shows the behavior of M and Z at low 
momenta. Arrows indicate the positions of the alleged 
pole in M, canceled by the zero of Z. 



XIV. OTHER PREDICTIONS 

The explicit model for the quark spectral function p{ui) 
allows for very simple and efficient evaluation of further 
interesting quantities. There is a whole bunch of pre- 
dictions, from which we only list a few. The results 
decouple into those involving the vector spectral func- 
tion, jOy(w), and the scalar spectral function, ps(w). As 
stresses throughout the paper, all predictions are made 
for the model working scale Qo, and the QCD evolution 
is needed if comparison to high-energy data is desired. 



Pion transition form factor 



evolution of the higher twist components. This point 
will be analyzed elsewhere. 



B. Pion light-cone wave function 



Next, we use Eq. H12.14JI in the expression for the pion 
light-cone wave function (|8.19|l . with the result 



3M3 



16tt(/4 + M2/4) 5 / 2 



6(x)6(l-x). (14.2) 



Passing to the impact-parameter space with the Fourier- 
Bessel transform yields 



*(x,6) 



27r / k±dk±®(x,k±)Jo(k±b) 
lo 



bMy 



exp 



M v b 



9(x)9(l-x). 
The expansion at small b yields 



(14.3) 



*M) = 1- 



M 2 b 2 
8 



M 3 b 3 
24 



x 9(x)9(l-x). 



(14.4) 



The average transverse momentum squared is equal to 



(*i> 



/ d 2 k 1 _ k ± ) _ d$(x, b) 

J d 2 k ± ^(x,k ± ) ~ ~ db 2 
M 2 



fc=0 



(14.5) 



which numerically gives (kj_) = (544 MeV) 2 (all at 
the model working scale Qo)- This value is not far 
from the result of the Nambu-Jona-Lasinio model with 
the Pauli-Villars regularization, which produces = 
(613 MeV) 2 [HI. The estimates from QCD sum rules 
yield smaller values: one gets (316MeV) 2 based on 
Ref. |H3, and (333 ± 40MeV) 2 based on Ref. 30]. One 
should note that these estimates are at Qo of the order of 
1 GeV, hence the QCD evolution is needed to compare 
to our model results [17|. 



According to the vector meson dominance model of the 
pion form factor, the transition form factor l|8.9[) becomes 



-^7T7*7*(Q 7^) 



2f v 1 
AN C Q 2 



log 



2M 2 



(1 + A)Q 



2 1 



2M 2 



M 2 



N C [4M* 



4Q 2 M 2 



(1 - A 2 )Q4] 



(1-A)Q 2 
(14.1 



As discussed in Sect. IVIIII this expression provides the 
twist expansion, with the result of Eq. (|8.17|) . but with 
VMD the log- moments of Eq. I|12.8|l . The analysis and 
comparison to the data of Ref. |52( requires the QCD 



C. Unintegrated pion structure function 



Due to relation 1)10.15(1 . exactly the same formulas as 
in the previous section hold for the unintegrated quark 
distribution function in the pion. Hence, we find 



q(x,k±_) = q(l - x,kx) 



(14.6) 



3M3 



167r(fc^ + M 2 /4) 5 / ; 



: 6(x)9{l-x), 



and the analogs of Eqs. i|14.3ll4.5fl . 
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D. Non-local quark condensate 



The applications of Sect. [XIV Bl and 1X1 V Ul involved 
the vector spectral density. The scalar spectral density 
(|13.1[) enters the evaluation of the nonlocal condensate 
of Eq. H2.18fl . For ds = 5/2 we find immediately the 
interesting and simple result 



Q(x) = exp 



M S V- 



(14.7) 



where x denotes the Minkowski coordinate. For higher 
values of ds the expression is multiplied by a polynomial 
in the x 2 variable. Note the nonanalyticity in the x 2 
variable in Eq. (|14.7[) as x 2 — ► 0. In the present model 
the moments of the quark condensate, Eq. (|2.20[) . are 
well defined for n < 2ds — 4, which means that for the 
preferred value of ds = 5/2 we may only consider the 
quark condensate itself, but not its moments. Neverthe- 
less, the coordinate representation of the non-local con- 
densate makes sense and can be given for the whole range 
of x 2 . 



E. Quark propagator in the coordinate 
representation 

One may also pass to the coordinate representation for 

introducing 



the A and B functions of Eq. (|2.7|) 

d A p 



B(x) 



(27T) 4 
d A P 



ip-x 



i>A{p) 



(2tt)* 

With straightforward algebra one finds 



ipx B(p) 



A(x) 



B(x) 



48 + 24My\/=x 5 - 6M 



2 Jl 



x e 



-M v 



-x 2 /2 
2 3-2d s (g q j 

N c T(d s - 2) 



967r 2 a 



(M s ^ 2 ~) ds - 2 K 2 ^ 



where x denotes the Minkowski coordinate. We have used 
dv — 5/2. In the limit of low x one recovers the result 
of the free theory, A{x) ~ l/(2n 2 x 4 ). For the preferred 
value of ds — 5 /2 we find that 



B(x) = (qq)/(4N c ) exp(-M sV /^2/2), (14.10) 
in agreement with Eq. l|14.7fl . 

XV. CONCLUSION AND FINAL REMARKS 

In the present work we have developed a chiral quark 
model, which tries to incorporate as many known fea- 
tures based on chiral symmetry and the partonic quark 



substructure of hadrons as possible. This approach, first 
unveiled in Ref. [l^, should be considered as a simple 
prototype of a construction which may be certainly im- 
proved in many respects. Taking into account the fact 
that the emerging picture is very encouraging, economic, 
and predictive, we believe that applications and exten- 
sions of the model deserve a thorough further investiga- 
tion. 

The key ingredient is the use of a generalized spec- 
tral representation for the quark propagator combined 
with the gauge technique for constructing the vertices in- 
volving quarks and currents. The generalized Lehmann 
representation implies analyticity properties of the quark 
propagator on the complex plane but positivity or real- 
ity are abandoned. In particular, the analytic continu- 
ation from the Euclidean to the Minkowski space back 
and forth becomes straightforward. We remind that this 
continuation is explicitly used when computing hadronic 
matrix elements in momentum space. The possibility 
of doing this continuation becomes extremely convenient 
when dealing with calculations of soft matrix elements of 
high-energy processes. In a purely Euclidean formulation 
one stays in the coordinate space and the extraction of 
the parton distribution functions or amplitudes is often 
limited in practice to the few lowest moments. 

The conditions used to constrain the spectral function 
are collected in Table I. Finiteness and factorization en- 
force the vanishing of the positive moments, while a num- 
ber of available experimental observables can be used to 
fix the values and signs of the negative moments and the 
log-moments. 

Instead of making a specific model for the propaga- 
tor based on extrapolations from the Euclidean region to 
the complex plane we devise an infinite set of spectral 
conditions based on the requirement that our model pro- 
duces finite hadronic observables. As a result the high 
energy behavior of certain matrix elements corresponds 
to a pure twist expansion, with no logarithmic behavior. 
This is consistent with the interpretation that the model 
is defined at a low renormalization point Qq. Thus, our 
d s (Ms\/ —x 2 /2), model results should be considered as initial condition for 
the QCD evolution, hence automatically incorporating 
(14. 9) the correct high energy radiative corrections to hadronic 
observables. The low renormalization point Qq is deter- 
mined by the analysis of two different processes: the pion 
transition form factor and the pion structure functions. 
We find compatible values in the range Qq ~ 300MeV, 
corresponding to a(Q Q ) / (2ir) ~ 0.3. For the model scale, 
the leading twist contribution for the pion structure func- 
tion and the pion distribution amplitude coincide and are 
equal to one, q(x) = ip^[x) = 1, regardless of the spec- 
tral function. At the same time the correct anomalous 
form factor for the process 7r° — > 77 is obtained. This 
resolves the conflict between proper low-energy and high- 
energy normalization of the pion transition form factor. 
The QCD evolution of both the PDF and PDA has been 
shown to provide a very reasonable agreement with the 
data analysis. 



(14.; 



r 2N3/2 
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TABLE I: Summary of conditions for the quark spectral function p(uS). 



Spectral condition 


Physical significance 


normalization 


po = 1 


proper normalization of the quark propagator 




preservation of anomalies 




proper normalization of the pion distribution amplitude 




proper normalization of the pion structure function 




reproduction of the large- N c quark-model values of the Gasser-Leutwyler coefficients 




relation M v = 24ir 2 f 2 /N c in the vector-meson dominance model 


positive moments 


Pi =0 


vanishing quark mass at asymptotic Euclidean momenta, M(Q 2 ) — > 


P2 =0 


finiteness of the pion decay constant, f-* 


pz = 


finiteness of the quark condensate, (qq) 


PA = 


finiteness of the vacuum energy density, B 


p„ = 0, n = 2,4,... 


factorization in the twist expansion of vector amplitudes 


/5 n = 0, n = 5, 7, . . . 


finiteness of (q(d 2 ) {n - 3)/2 q) 




factorization in the twist expansion of the scalar pion form factor 


negative moments 


p-2 > 


positive value of the quark wave-function normalization at vanishing momentum, Z(0) > 


P-l/P-2 > 


positive value of the quark mass at vanishing momentum, M(0) > 


P-n 


low-momentum expansion of correlators 


positive log-moments 


P2 < 


f 2 = -N c /(4n 2 )p' 2 


P3>0 


negative value of the quark condensate, (qq) = —N c /(4ir 2 )p' 3 


P4 > 


negative value of the vacuum energy density, B — ~N c /(Att 2 )p' 4: 


P5 <0 


positive value of the squared vacuum virtuality of the quark, \ 2 — —p'zlp'z 


Pn 


high-momentum (twist) expansion of correlators 



The pion form factor does depend on the spectral func- 
tion. Further interesting analytic relations can be ob- 
tained by determining the even contribution to the spec- 
tral function from the requirement that the pion form 
factor has the vector meson dominance form, which is 
known to describe experimental data in the available mo- 
mentum range very well. As a result, the vector meson 
mass becomes proportional to the pion weak decay con- 
stant. The vector contribution to the spectral function 
exhibits a pole at the origin and non-integrable branch 
points at plus or minus the vector meson mass. As a 
consequence the spectral function must be interpreted as 
a distribution centered at the branch points. 

The form of the spectral function suggested by the vec- 
tor channel can be also used with minor modifications for 
the scalar channel. This way the full quark propagator 
can be studied and its analytic structure analyzed: there 
are no poles all over the complex plane, only cuts located 
at the branch points of the spectral function. Moreover, 
the asymptotic behavior in the Euclidean region reflects 
this cut structure by a half-integer power fall-off of the 
quark mass function in the squared Euclidean momen- 
tum Q 2 . The recent lattice data 0, 0] set this half 
integer power to 3/2 quite unambiguously. Fixing this 



power leaves only the quark mass function at the origin 
(the constituent quark mass) and the quark condensate 
as free parameters. The fit to the data for the mass func- 
tion, M(p 2 ), is very good and the values for both the 
condensate and the constituent quark mass at p 2 = 
agree with other estimates. The fit to the quark wave- 
function normalization, Z(p 2 ), is not nearly as good as 
for M(p 2 ), leaving room for improvement. 

Possible extensions of the general model involve the in- 
clusion of more general solutions of the Ward-Takahashi 
identities. The meson-dominance model can be improved 
by studying more general forms of the scalar quark spec- 
tral function. Moreover, the inclusion of finite quark 
masses would allow the extension of the present model 
to the complete pseudoscalar octet. These issues are un- 
der investigation. 
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APPENDIX A: ONE-LOOP INTEGRALS 

1. Two-point integral 

The two-point one loop integral regularized in 4 + e 
dimensions is 



i r d*k 



i 



i J (2tt) 4 k 2 - uj 2 + i0+ (q - k) 2 -lu 2 + M+ 

£-1 



4w 2 



1 



16tt 2 



1 log ■ 

1 ./I 



1 



with 



/(0 '" ) = -ii(7 + log( ^ 2) 



+ I(0,w), 
(Al) 

(A2) 



We also introduce 

I(q 2 ,u) = I(q 2 ^)-I(Q,Lo). (A3) 
In the Feynman parametric form we equivalently have 
1 



dup(uj) 



x)q- 



(A4) 



(4tt) 2 

x / dxlog [lo 2 + x(l 
Jo 

The imaginary part yields 

IIm/( 9 ^.) = I ^L^l~Z^-4 W 2 ). (A5) 
Thus the once-subtracted dispersion relation, 



2 roo 



dt Im/(i, to) 
,2 T t - q 2 - iQ + ' 



(A6) 



holds. The asymptotic behavior for large Euclidean — q 2 
is 



1 



1 



16tt : 



{2 log(- 5 7^ 



[log(-,> 2 ) + 1] + ^ 
At low q 2 we have 



log(V/^ 2 ) 



Hi 2 ,") 



i 



q" 



16tt 2 6uj 2 60cj 4 420w 6 



(A7) 
...}. 

(A8) 



2. Three-point integral 

The three-point one-loop integral is defined as 

d 4 k 1 



K{{qi - q 2 ) 2 ,ql,ql,uj) = ~ / 

i J (2tt 



) 4 k 2 ~LU 2 + i0+ 



1 



1 



(k - qi) 2 -uj 2 + i0+ (fc - q 2 ) 2 -lu 2 + i0+ ' 



(A9) 



We analyze it with the dimensional regularization, and 
for the case where the virtuality of one of the external line 
vanishes, (qi — q^) 2 — (massless pion). We immediately 
find the result 



JsT(0,0,0,w) 



1 1 



16tt 2 2lo 2 

The following Feynman parameterization is useful 
1 



(A10) 



- = 2/ dxi dy 



l [xya + x(l — y)b + (1 — x)c] 3 

(All) 

Carrying the momentum integration and introducing z = 
2y — 1 yields 



K(0, 



1 + A 2 1 - A 



-q , 



2 * ' 2 

X 



q 2 M 



dx dz x 



(1 - x)x(l + Az)q 2 - 2uj 2 ' 
At A = we find 



KiO^q 2 ,^ 2 ^) 
1 2 



(A12) 



(A13) 



dog 



1 + y/1 - 8w 2 /q 2 



16tt 2 q 2 ^/l-8LU 2 /q 2 1 - y/1 - 8w 2 /q 2 
At low q 2 the expansion is 



(A14) 



1 



16?r 2 
and an high g 2 



A'(0,^) 



1 



2uj 2 24uj 4 24Qu 4 



16tt 2 



21og(-2^ 2 /9 2 )- 



-8 [log(-2c 2 /g 2 )+l]^ + . 
The integral over x in Eq. (|AT2)l yields 



1 + A , 1 - A 



-9 > 



-9 = 



1 



log 



1+8 



K(0, 

16-7? ( I f:l:kr/" n l -.s' 
s= v / l-8o; 2 /( (? 2 (l + Az)). 



(A15) 



(A16) 
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At low q 2 we have 



We choose the reference frame of the target at rest 



(A17) p = ( m , ± , 0), q = (q ,0 ± ,q 3 ), q 2 = -Q 2 = q 2 - q 2 . 



16tt 2 



q 2 (A 2 + 3)g 4 



1 



2oj 2 24w 4 720u 4 



The large-g 2 expansion produces 



K(0, 



^4 2 ,^4» (A18) 

1 f 1 , I 1 ] ( 2co 2 
16tt 2 J_ x \q 2 (l + Az) ° S \ q 2 (l + Az) 



Auj 2 



q\l + Azf 



log 



2u 2 



q 2 (l + Az) 



1 



The integral over z in Eq. I|A16(I gives finally the simple 
general result 



K(0. 



l+A o I- A 



-q i 



( f) 



l 



327rMg 2 



log 



log 



1 + 5- 
1 - S_ 



(A19) 

2" 



S± = Vl-8 W 2 /(<z 2 (l±A)). 



(B3) 



One gets then 



<7o 



<7:S 



Q 2 



Q 2 I | 4m^ 2 _Q 2 



(B4) 



7^2 ' o + ™"* x + 



2m 7T x 

In the light-cone coordinates, defined as 

k — k ~\~ k . k — k ~~ k . — ^/c . k 1^ 



one obtains 



(B5) 



«?+ = q° + q 3 ^-^- 



q = q° - q 3 -> -m„X, 



(B6) 



APPENDIX B: DISCONTINUITY IN THE 
BJORKEN LIMIT 

Let us consider the one-loop function 



T(p,g) = i 



d A k 



1 



(2tt) 4 [(* - p) 2 -lu 2 + i0+} 2 X ^ 



and also 

6+ [(k-p-q) 

5+ [k 



2 u?\ 



~^2~ s [k - (1 - a;)m 7r ] , 



'~ 2 ^ 2 1 



1 



771^(1 — £k) 



171^(1 — x) 



(B7) 



fc 2 - lu 2 + i0+ (q - k - p) 2 - uj 2 + i0+ ' 

The discontinuity in the s = (p + q) 2 channel may be 
computed through the Cutkosky rules, 



DiscT M = */^ 



1 



r) 4 [{k-pf -lu 2 + i0+} 2 
{-2iri) 2 5+{k 2 -uj 2 )8+[{q-k-p) 2 -uj 2 ]. (B2) 



Thus finally we get 



ImT(p, q) 



(B8) 



x(l — x) f d 2 k 
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[67] Our conventions are as follows: We take 



for 



going and e w ' x for out-going particles. For free massless 
particles the irreducible functions are normalized accord- 
ing to Vy a = 7 M %- and F^ a = 7 M 7 5 % L . The convention 
for the Dirac matrices and the metric tensor is the one 
of Ref. 

The use of the dimensional regularization guarantees 
gauge invariance, with no further subtractions, hence for 
simplicity we use it throughout the paper in the applica- 
tions with vector and axial currents. 
For instance, a modification of the vector vertex by pro- 
viding additional tensor coupling is capable of curing the 
problem of the second Weinberg sum rule. This impor- 
tant issue will be studied elsewhere 

If one extrapolates pv(u>) to the real axis, as suggested 
by the standard Lehmann representation, Eq. (12. 211 . there 
appear, for dy > 1, end-point non-integrable singulari- 
ties at the branch points uj = ±My/2. If one insists on 
a real spectral function, one may do so but then one has 
to proceed by analytic continuation in dv, or derivation 
with respect to My after computing the oj-integral. Al- 
ternatively, this is equivalent to a distributional interpre- 
tation of the spectral function p{uj) and its (generalized) 
derivatives using the well known distribution x% (see e.g. 
the classical work |66jl for a rigorous discussion). 
Obviously, M(p 2 ) and Z(p 2 ) are gauge-dependent quanti- 
ties. The data of Ref. [181 Hfl| correspond to the Landau 
and Laplacian gauges. A natural question arises where 
our results incorporate the choice of the QCD gauge at 
the microscopic level. We notice two sources of ambigui- 
ties in our approach: the arbitrariness of the transverse 
terms in the solutions to the Ward-Takahashi identities, 
as well as the choice of the scalar quark spectral func- 
tion. It remains to be seen how these choices reflect the 
underlying QCD gauge. 



